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ABSTRACT 



> 
On 

■^j- ■ We discuss if-band (1.65 //m) near-infrared photometry of the central 9 h 

Mpc 2 of Abell 1644 to a limiting M H ~ MJj+3 (throughout this paper H = 100/t 
km s _1 Mpc -1 ). There are 861 galaxies in the photometric survey region. We also 
measured radial velocities of 155 galaxies; 141 of these are cluster members within 
2.44 h~ l Mpc of the cluster center. The completeness limit of the spectroscopic 
survey is H ~ 13 (M# ~ M^). The cluster velocity dispersion of a ~ 1000 
km s _1 remains constant out to the limiting radius. We find no evidence for 

d ■ 

substructure in the cluster. The cluster mass within R = 2.4 h' 1 Mpc is 7.6 ± 



1.3 x 10 14 h- 1 M Q . 

We compute the cluster luminosity function; the Schechter parameters a 



X 
& . 

-1.14 ± 0.08 and M* H = -24.3 ± 0.2 (with h = 0.5) agree well with other H- 
band luminosity functions. From the virial theorem and the caustic method 
we compute one of the first mass-to-light ratios at H; the result is M/Lh = 
82-127/iM /L within 1.5/r 1 Mpc. This ratio corresponds to 374-579/iM /L 
at R. The agreement of our IR measurement with previous M/L determinations 
indicates that at low redshift dust and young stellar populations may produce 
only negligible systematic errors in optical mass-to-light ratios. 

Subject headings: galaxies: clusters: individual (A1644) — galaxies: luminosity 
function, mass function 
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1. INTRODUCTION 

Observations of galaxy clusters are a useful tool for constraining models of structure 
formation and of determining cosmological parameters. Cluster mass-to-light ratios, tradi- 
tionally measured at optical wavelengths, are important in estimating Q (c.f. Bahcall 1977; 
Bahcall, Lubin, & Dorman 1995; Carlberg et al. 1996; and references therein). The lumi- 
nosity functions (LFs) of clusters can also be compared to field LFs to investigate the effects 
of galaxy evolution in the cluster environment. 

Infrared observations of clusters have several advantages over optical measurements. IR 
emission is less affected than visible light by absorption and scattering due to dust, which 
may be interstellar or intracluster (Zwicky 1962; Boyle, Fong, & Shanks 1988) in origin. 
Galactic infrared emission is also dominated by the older stellar population and is relatively 
insensitive to short-lived blue stars and changes in the star formation rate (Bruzual & Chariot 
1993). Infrared radiation should therefore be a good tracer of the total stellar mass of the 
cluster. 

Until now, the small fields of view of infrared CCD arrays have impeded deep large-area 
near-infrared photometry of clusters. For example, Andreon & Pello (2000) measure the 
H-bsaid LF of the Coma cluster to a faint limiting magnitude (M^ + 6), but in a projected 
area of 2% of the size of our survey. The 2MASS project surveys the entire sky but is almost 
2 mag shallower than our data at H. Cluster surveys at K are slightly more common (De 
Propris et al. 1999; Trentham & Mobasher 1998; Barger et al. 1998), but most suffer from 
similar limitations. 

The iZ-band data we discuss here, covering a projected area of 9 h~ 2 Mpc 2 and reaching 
Mh ~ Mh + 3, make Abell 1644 (A1644) one of the best-studied clusters at near-infrared 
wavelengths. We constrain the IR LF and compare its shape and parameters to those 
for other cluster and field surveys. The addition of extensive spectroscopy enables us to 
comment on several kinematic issues. We examine the properties of emission and absorption 
galaxies, test for substructure, and compute the infrared mass-to-light ratio of the cluster. 
This measurement, Mj Lh ~ 82 — 127 h M Q /L Q , is consistent with observed optical values 
and is one of the only existing cluster M/L determinations in the infrared. 

In Section 2 we describe our photometric and spectroscopic observations. Section 3 
contains an analysis of the kinematic properties of the cluster. In Section 4 we compute the 
LF and mass-to-light ratio, and we summarize in Section 5. 
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2. OBSERVATIONS 

Abell 1644 is a nearby rich cluster with a cD galaxy (a = 12 h 57 m ll s .6, 5 = -17°24'34".4 
J2000, cz = 14,233 ± 33 km s -1 ) at its center. Our H-b&nd photometric observations cover 
a 3 x 3 hr 2 Mpc 2 region surrounding the cluster core to limiting H ~ 16. We also measured 
redshifts for 155 galaxies within 2.44 /i -1 Mpc of the cluster center. The combined sample 
of 123 cluster members with both spectral and photometric data is complete to H ~ 13. 

2.1. Photometry 

We obtained i/-band photometry covering a 36' x 36' region centered on a = 12 h 57 m 35 s 
5 = -17°27'21" on 1997 April 16-19. We used the CIRIM infrared camera on the 1.5 m 
telescope at Cerro Tololo Inter- American Observatory (CTIO). The scale for the 256 x 256 
array is 1".137 pixel -1 , and the duration of each exposure was 6 x 10 s . We operated the 
instrument in a mapping mode which acquired pairs of dithered images offset by 10". We 
then covered the survey region with a tiling of dithered pairs which overlap by 40". In 
addition, we surveyed two off-cluster fields of sizes 0.22 deg 2 and 0.04 deg 2 to ascertain 
background statistics. 

On the nights of 1999 April 25-29 and 1999 May 26 we used the Stelircam IR camera 
on the 1.2 m telescope at Fred L. Whipple Observatory (FLWO) on Mt. Hopkins, Arizona 
to expand our infrared survey. We operated the 256 x 256 instrument (1".2 pixel -1 ) in a 
mapping mode similar to that used at CTIO. Of the six nights at FLWO, two were not 
photometric. On the remaining nights we gathered data over a 0.99 deg 2 region bordering 
the CTIO observations (bringing our total coverage in the cluster to 1.35 deg 2 ), and we 
increased the sizes of our off-fields to 0.41 deg 2 and 0.44 deg 2 . We also imaged a small 
area in the center of the cluster, overlapping the CTIO data, to assess systematic differences 
between the two surveys. 

For each photometric night, we first used NO AO IRAF 2 commands to flat-field and 
sky-subtract each frame. We calculated the zero point of each night's magnitude scale by 
fitting a linear solution to several standard star images and correcting each object frame for 
airmass. We then used the SExtractor software package for source detection and photom- 
etry. In addition to the usual aperture and isophotal magnitudes, SExtractor calculates a 
composite MAG_BEST. MAG_BEST is useful in crowded fields and is equal to the Kron 
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magnitude (Kron 1980) if no neighbor can bias the Kron aperture magnitude by more than 
10%. Otherwise it is the corrected isophotal magnitude (Bertin & Arnouts 1996). We adopt 
MAG_BEST for our subsequent data analysis. To demonstrate the consistency of the pho- 
tometry, we plot the magnitude differences as a function of MAG_BEST for objects observed 
in successive dithered frames on 1997 April 17 (Figure 1). The rms error in the photometry 
is Am rms ~ 0.12 mag at H = 16. 

The SExtractor software also includes an algorithm to separate stars and galaxies based 
on morphology. Each object identified has a value CLASS_STAR which ranges from zero 
for galaxies to one for stars. Ambiguous objects lie somewhere in between. We examine 
the images manually and find that CLASS_STAR accurately separates stars and galaxies 
brighter than H ~ 14.5; fainter objects are indistinguishable by eye in these images. 

An examination of the CTIO and FLWO data in the overlap region reveals that the 
seeing is better in the CTIO images. However, there is little systematic difference in the 
photometry. The average magnitude difference between objects observed at FLWO and 
CTIO is Am = 0.04 mag, and the rms scatter is Am rms = 0.15 mag. The SExtractor 
star/galaxy separation gives consistent results for objects observed at both sites. In the 
overlap region, 93% of the objects with H < 14.5 identified as galaxies in the FLWO data 
are also classified as galaxies in the CTIO data. 

Table 1 shows photometric data for the 861 galaxies in our infrared survey. For H < 13 
the galaxy spectroscopic survey is complete; for 13 < H < 16 we list all objects with 
CLASS_STAR < 0.5. Columns (1) and (2) give the right ascension and declination of each 
galaxy; column (3), the apparent H magnitude (MAG_BEST); and column (4), the redshift 
if measured. 



2.2. Spectroscopy 

We acquired optical spectra for 155 galaxies in the region of Abell 1644 during 2000 
January- June with FAST, a high throughput slit spectrometer mounted on the 1.5 m tele- 
scope at FLWO. The FAST setup includes a 300 groove mm" 1 grating and 3" slit, and yields 
spectra which cover 3600-7600A at a resolution of ~ 6A. Although Dressier & Shectman 
(1988a, hereafter DS) measured 103 redshifts in the region of Abell 1644, we re- acquired 
spectra of 91 of their objects to assure uniformity throughout our data set and to assess the 
spectral type of the galaxies. The remaining twelve DS galaxies are either too faint for FAST 
redshift measurements or are outside our survey region. We selected 64 additional targets 
on the basis of the DS catalog, digitized sky survey images, and our own infrared survey. 
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We obtained a reduced spectrum for each object using IRAF. We made bias, flat- 
field, and illumination corrections to each frame and then used the IRAF XCSAO task to 
obtain redshifts via cross-correlation with standard emission and absorption templates. The 
resulting spectroscopic sample is complete to H ~ 13 in the region where we have i7-band 
photometry. We tested the differences between the FAST and DS observations for the 91 
common galaxies. The FAST velocities are systematically lower by ~ 39 km s _1 , and the 
rms velocity difference is Av rms = 80 km s _1 . These offsets are consistent with the estimated 
velocity errors of roughly 30 km s" 1 (FAST) and 45 km s -1 (DS). 

Table 2 lists properties of the 155 galaxies with measured redshifts. Columns (1) and 
(2) give the right ascension and declination; column (3), the radial velocity; column (4), 
the error on the radial velocity; column (5), the spectral type (either emission or absorption 
as discussed in Section 3.1); column (6), MAG_BEST (if available); and column (7), the 
reference number from the DS catalog for galaxies in the DS list. 

2.3. Cluster Sample 

Figure 2a shows the measured velocity distribution. A1644 is remarkably isolated in 
redshift space. A few background galaxies are visible at redshifts up to cz ~ 65, 000 km 
s _1 . There is also evidence of foreground structure in the range 2500 < cz < 10,000 km s _1 . 
Figure 2b shows the 141 cluster members with redshifts between 10,000 and 20,000 km s _1 . 

The mean radial velocity of the cluster is cz = 14, 295 ± 93 km s _1 with a dispersion 
of a = 1108jZgf km s -1 , computed at the 68% confidence limit as in Danese, De Zotti, & 
di Tullio (1980). The radial velocity of the cD (14,233 ± 33 km s" 1 ) is coincident with the 
overall cluster mean. The cluster sample is isolated from foreground and background galaxies 
by velocity gaps of 3244 km s" 1 on the low end and 4589 km s" 1 on the high end. The sizes 
of these gaps are 2.9a and 4.1(7, respectively. There are 123 cluster members within the 
region of infrared photometry; Figure 2c shows the velocity distribution of these galaxies. 

3. KINEMATICS 

3.1. Emission and Absorption Galaxies 

We separate the galaxies into emission (Em) and absorption (Abs) systems. Emission 
galaxies are those objects with a larger R value when cross-correlated with the emission 
template; all other galaxies are classified as absorption. Figure 3 shows the positions of the 
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141 cluster members according to spectral type. There are 18 emission and 123 absorp- 
tion systems. The absorption galaxies are more centrally condensed; the median projected 
distances of the Abs and Em samples from the cD galaxy are 20'. 1 and 27'. 5, respectively. 
This difference is only marginally statistically significant; a rank-sum test applied to the two 
samples gives a 97% confidence that the Em and Abs systems do not come from the same 
parent spatial distribution. 

Figure 4 shows the radial velocity distributions of the Abs and Em samples. The 
absorption galaxies have cz = 14, 363 ± 98 km s" 1 and a = 1095^ km s -1 . The dispersion 
of the Em sample (a = 11231^ km s _1 ) is comparable, but the mean velocity of cz = 
13, 800 ± 273 km s -1 is lower. However, a rank-sum test gives only a 96% confidence that 
the Em and Abs galaxies follow different velocity distributions. 



3.2. Substructure 

Dressier & Shectman (1988b) define a test parameter 5 which measures the deviation of 
a subgroup's local mean velocity and dispersion, Vi oca i and <Ji oca i, from that of the cluster's 
overall values (v and a): 

5 2 = —^[{Vlocol ~ vf + ((Tlocal - <j) 2 ], (1) 

where n is the number of galaxies in the subgroup. The cumulative deviation is A = YliLi <^> 
where N is the number of cluster members and for each galaxy we calculate Si using the 
galaxy and its n—1 nearest neighbors. For a Gaussian velocity distribution with only random 
fluctuations, A ~ iV; the presence of substructure can cause A to be significantly higher 
than N. 

Dressier & Shectman (1988b) analyze 92 members of A1644 with n — 11 and find 
evidence for substructure in the southeastern quadrant of the cluster. To test the significance 
of this structure, Dressier and Shectman randomly shuffle the velocities of the galaxies to 
create roughly 1000 model clusters. They calculate the value of A for each model and find 
that only 2.7% of the Monte Carlo clusters have A > A observe d. 

We duplicate the Dressier and Shectman analysis for our larger sample of 141 cluster 
members. Allowing n to range from five to 105, we calculate the A statistic for each subgroup 
size. We also create 5000 Monte Carlo models for each value of n and compute the probability 
that A > A observed . Figure 5 shows the results. 

With n = 11, we find that A is greater than or equal to A obS crvcd in 17.5% of the 
simulated clusters. Thus our larger sample, when compared to the Dressier and Shectman 
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catalog, no longer provides evidence for substructure on this scale. The best evidence for 
substructure occurs with n = 6 and n = 67. However, these detections are not statistically 
significant, with probabilities of A > A observed of 4.3% and 5.0%, respectively. 



Figure 6 shows an azimuthally averaged radial profile of the cluster velocity distribution 
for both Abs and Em galaxies. We compute the velocity dispersion of the cluster as a 
function of distance from the cD galaxy and plot the results in Figure 7. We choose the 
cD as the geometric center of the cluster because of the lack of evidence for substructure. 
Moving outward from the cD, we calculate the velocity dispersion of successive groups of 11 
galaxies arranged in order of increasing radius. The error bars in the plot represent the 68% 
confidence level for uncorrelated points. 

The velocity dispersion of the cluster remains roughly constant for R < 1.9/i -1 Mpc; 
at this radius it nearly doubles. The increase at large radius is extremely sensitive to the 
presence of three high- velocity galaxies circled in Figure 6. When we exclude these galaxies 
(Figure 7, bottom plot) the velocity dispersion at R ~ 1.9 Yr 1 Mpc decreases by a factor of 
three. We suggest that the three high- velocity galaxies may be outliers and that a is roughly 
constant as a function of radius. Treating these objects as background is consistent with 
the results of our caustic analysis (Section 3.4), and also with Dressier & Shectman (1988b), 
who select cz < 18, 000 km s _1 as a criterion for membership in A1644. 



We calculate the mass of A1644 in two ways. Our first mass estimate depends on the 
"caustic" technique of Diaferio & Geller (1997). A hierarchical clustering model predicts 
the existence of two caustic curves with amplitude A(R) approximately equal to the escape 
velocity from the cluster at radius R. Diaferio (1999) shows that A(R) is related to the mass 
of the cluster interior to R: 



We use the techniques of Diaferio (1999) with smoothing parameter q = 25 to calculate 
the caustics of A1644 (Figure 8). Of the 141 presumed cluster members, 127 lie within the 
caustics; the fourteen remaining galaxies may in fact be outliers. Equation 2 automatically 
excludes these outliers from the caustic mass determination. 



3.3. Velocity Dispersion Profile 



3.4. 



Cluster Mass 




(2) 
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Our second mass estimate uses a virial estimator (Binney & Tremaine 1987): 



3ttN Eti^p. 
2G T N V \R - R 



M est = ~^T ^ N JT l= " „, (3) 



where t> Pi j is the radial velocity of each galaxy with respect to the cluster mean and Ri is the 
galaxy's position relative to the cD. This estimator assumes that the galaxies are embedded 
in a diffuse dark matter distribution, and that the spatial arrangement of galaxies traces the 
dark matter. 

The virial mass is very sensitive to outliers. In Figure 9 we first plot M V i r i a i using the 
141 galaxies with 10,000 km s _1 < cz < 20, 000 km s _1 . The virial mass decreases by a factor 
of two or more, even at large radius, when we include only those galaxies which lie inside 
the caustics. For comparison, we plot M caustic for three different values of the smoothing 
parameter; the caustic mass is much more robust. 

Equation 3 overestimates the true value of M because we have not subtracted the virial 
surface term (The & White 1986). When the surface term C(r) is included, the corrected 
virial mass is 



M cv (< b) = M est {l - C(b)} = M e Jl-4nb 3 f-^iU 2 }, (4) 



where p(r) is the radial density distribution, ay is the radial component of the velocity 
dispersion, and er(< b) is the integrated velocity dispersion within the limiting radius b 
(Girardi et al. 1998). 

Calculation of the surface term requires knowledge of both the velocity dispersion profile 
and the density distribution. While a(r) can be constrained from the data, the density profile 
is unknown and is, in fact, the very quantity we are attempting to correct via equation 4. 
To approximate p(r) we fit a Navarro, Frenk, & White (1995, 1996, 1997, hereafter NFW) 
mass profile to M caustic . The NFW functional forms of the mass and density profiles are 

M(< r) = 47i5 c p c rj [log(l + r/r c ) - r J r ° 1 (5) 

1 + f/fc- 



P(r) = TZTZXlf^r^' (6) 



and 

8cPc 

(r/r c )(l + r/r c f 

where 8 C and r c are model parameters and p c is the critical density of the universe. We fit 
equation 5 to the caustic mass profile and substitute the best-fit parameters into equation 6 
to recover the density distribution. We then use the fitted p(r) to estimate the virial surface 
term via equation 4. 
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To estimate C(r) one must also assume a value of the anisotropy parameter j3{r) = 
1 — oj/of • Girardi et al. (1998) show that clusters with a flat velocity dispersion profile are 
best described by models with f3{r) = 0. We adopt this value of j3{r) and compute the virial 
surface term. The overall trend of the surface term is decreasing with radius. However, 
C(r) fluctuates on small scales because of its dependence on <r(r). At r = 1.5 hr x Mpc, 
C(r) -0.2 ±0.1. 

We use a statistical "jackknife" procedure (Diaconis & Efron 1983) to calculate errors 
in the virial mass, AM/M ~ 2 — 4% for R = 0.5 — 1.5 h^ 1 Mpc. Systematic sources of error 
in the virial mass dwarf this formal jackknife estimate; the surface term correction is highly 
uncertain due to the assumptions regarding p(r) and (3{r). In addition, the virial estimate 
depends strongly on the particular galaxy sample. 

Figure 10 shows our final estimates of the caustic and virial mass profiles. For consis- 
tency, we base both estimates on the same set of 127 galaxies enclosed by the caustics. Errors 
in M caust ic represent uncertainty in the location of the caustics. The virial mass incorporates 
the surface term correction. We show M viria i as an envelope whose maximum and minimum 
values signify the range of masses consistent with the jackknife errors and the uncertainty 
in C(r). With our adopted errors, the caustic and virial mass estimators are consistent at 
the 2a level at all radii. 

Figure 10 also includes two X-ray mass determinations from Ettori, Fabian, & White 
(1997). They arrive at mass estimates for A1644 based on Einstein Observatory IPC X-ray 
data. Their deprojection analysis yields M dpr = 1.05 x 10 14 h~ x M Q within 0.305 h~ l Mpc. 
They then extrapolate the deprojection mass to R500, the radius where p/p cr u = 500; the 
extrapolated mass is M500 = 4.06 x 10 14 h^ 1 M & within 0.85 h~ l Mpc. Our caustic mass 
estimate is consistent with the X-ray data at both radii. 



4. PHOTOMETRIC PROPERTIES 

4.1. Luminosity Function 

Figure 11 shows the differential iif-band luminosity function of A1644. In the bins with 
complete spectroscopy (H < 13) redshifts determine cluster membership. At the faint end 
we estimate the number of field galaxies because we do not have spectra. We model the 
number of field galaxies in each magnitude bin by a power law of the form 

N(rm <m<m 2 ) = C / 10 a67m rfm, (7) 

J mi 
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where C is a normalization constant. We assume the background slope of 0.67 determined 
by Gardner, Cowie, & Wainscoat (1993) from a compilation of ZT-band field surveys. As our 
best estimate of the background, we normalize the power law by integrating equation 7 from 
H = — oo to H — 14 and equating the left hand side with the number of field galaxies in the 
direction of the cluster; there are 5 ± 2 field galaxies deg -2 with H < 14 in our photometric 
region. Our off-cluster survey contains 39 ±6 bright galaxies per square degree, and Gardner 
et al. (1993) and Szokoly et al. (1998) observe roughly 13 ± 3 and 23 ± 6 galaxies deg -2 
in this magnitude range. Although these field counts are statistically consistent with our 
estimate, the large variations in the normalization can produce a rising, flat, or declining 
faint end of the A1644 cluster LF. 

After subtracting the background fit normalized by our redshift survey from the total 
number of faint galaxies identified by the SExtractor CLASS_STAR parameter, we charac- 
terize the data with a Schechter luminosity function. We fit to the Schechter (1976) form of 
the luminosity function 

n e {L)dL = n*(L/L*) a exp(-L/L*)d(L/L*), (8) 

where a is the faint-end slope, L* is a characteristic luminosity, and n* is a normalization 
constant. In terms of absolute magnitude M the luminosity function is 

N e (M)dM = kn* exp{[-A;(a + 1)(M - M*)] - exp[-k(M - M*)]}dM, (9) 

where k = (In 10/2.5) and M* is the absolute magnitude corresponding to L* (Kashikawa et 
al. 1995). At the mean redshift of our cluster M H = m H — 35.78 + 51og/i. 

We fit a Schechter function to our data by minimizing the quantity 

X , B E W) -N.m? t (10) 

where N(M,j) and iV e (Mj) are the observed and fitted numbers of galaxies in the ith mag- 
nitude bin and Oi is the variance in the ith measurement. With respect to the maximum 
likelihood method (Sandage, Tammann, & Yahil 1979), also commonly used to fit luminosity 
functions, a x 2 fit has the disadvantage that the data must be binned, but it is advantageous 
because it provides a measure of the goodness of the fit. We assume a Poissonian variance 
in the total number of galaxies and in the number of background galaxies in each bin. 

We determine the best-fit Schechter function (solid curve in Figure 11); the best-fit 
parameters are a = -1.14 ± 0.08, M* = -24.3 ± 0.2, and n* = 114 ± 22 (here and in 
subsequent discussions of absolute magnitudes we take h = 0.5). The Schechter function is a 
good fit to the faint end of the LF and slightly underpredicts the number of bright galaxies. 
The minimum value of x 2 is 6.3 with 8 degrees of freedom. 
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Our best-fit values of a and Mjj are in good agreement with a survey of A1644 in the 
visible (Kashikawa et al. 1995). They fit a Schechter LF (with a fixed at -1.25) to the R- 
band luminosity distribution and find Mjt, = —21.50. With an expected elliptical color of 
(R — H) ~ 2.5, the characteristic magnitude at H is ~ —24.0, which agrees very well 
with our measurement. 

There are several A"-band field luminosity functions. Parameters derived from these 
surveys span roughly —1.0 < a < —1.3 and —25.1 < M* K < —23.6 (Mobasher, Sharpies, & 
Ellis 1993; Szokoly et al. 1998; Loveday 2000). Our values for A1644 are consistent with this 
range, assuming H — K = 0.25 as above. Parameters derived from iT-band studies of clusters 
are similar. De Propris et al. (1999) find M* K ~ —24.9 ± 0.5 with a fixed at -0.9 for clusters 
in their z = 0.15 bin. A separate survey of five X-ray luminous clusters finds a composite 
a = —1.38 ± 0.24 (Trentham & Mobasher 1998). Our results are consistent with these low 
z counts for typical early-type galaxy color, but not with Barger et al. (1998), who find an 
average Mj*- = —25.4 ±0.1 for 10 clusters with 0.31 < z < 0.56. The discrepancy could be 
due to luminosity evolution, although Barger et al. claim that their M* K is not significantly 
brighter than the z ~ values. 

At H fewer data are available. An #-band study of a 0.1 deg 2 region of the Coma 
cluster finds a = -1.3 ± 0.2 and M* H = -24.6 ± 1.0, with a dip in the LF at M H = -22.2 
(H ~ 13.5) (Andreon & Pello 2000). We estimate the errors on their values of a and by 
examining the 68% confidence contours in their Figure 3. An imaging study of A548 over a 
much larger area (0.9 deg 2 ) reveals a similar dip at H ~ 13 — 13.5 (S. J. Kenyon et al., in 
preparation). The number of galaxies in our survey is too small to evaluate the presence of 
a dip at H ~ 13.5. The counts in the 13 < H < 13.5 and 13.5 < H < 14 bins (55 ± 7 and 
68 ± 8, respectively) are very similar; the LF in this region is thus consistent with either a 
flat, rising, or declining slope. 



4.2. Mass-to-Light Ratio 

We determine the total f/"-band luminosity of A1644 by extrapolating the luminosity 
function. For a magnitude-limited survey described by a Schechter function the observed 
fraction of the total light is given by T(a + 2, L min /L*)/T(a + 2), where T(x,y) is the 
incomplete Gamma function and L min and L* are the luminosities corresponding to the 
completeness limit and M*. With a completeness limit of H = 16, the fitted Schechter 
parameters imply that we observe 90.6% of the total light. Our photometric survey covers 
a square field 9 h~ 2 Mpc 2 in area. The photometry is therefore complete within 0.5(9 h~ 2 
Mpc 2 ) 1 / 2 = 1.5 h' 1 Mpc of the center of the field. Inside this radius the total luminosity 
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is L# = 7.6 ± 0.4 x 10 12 h 2 L Q ; the quoted error in L H reflects the uncertainty in the LF 
parameters. 

For comparison, we integrate the visible LF in a smaller 0.82 deg 2 region of A1644 
(Kashikawa et al. 1995). Using their values of a, M R , and M^ m , the total i?-band light 
within 1.1 hr x Mpc is L R ~ 1.1 x 10 12 hr 2 L & . Within the same radius we compute the H- 
band luminosity and use typical galactic colors and (R—H) & = 0.85 to convert H photometry 
to R. The result is L H = 5.1±0.3x lO 12 /i~ 2 L which implies L R ~ 1.12±0.06x lO 12 /i~ 2 L , in 
excellent agreement with Kashikawa et al. The agreement implies that the optical luminosity 
of A1644 is only slightly affected by factors such as dust extinction and recent star formation 
activity. 

Figure 12 shows the radial M/Lh profile. The mass-to-light ratio is flat for R > 0.4 h~ x 
Mpc, suggesting that the dark matter fraction is constant outside this radius. The increased 
M/L at R < 0.35 h~ x Mpc is a result of a depressed luminosity near the cluster core. At 
the limiting photometric radius of 1.5 h~ x Mpc the caustic and virial mass-to-light ratios are 
M/L H = 82 ± 13 h M & / L Q and M/L H = 127 ± 26 h M Q /L Q . 

For A1644, Girardi et al. (2000) obtain M/L B] ~ 250 h M & /L & within 1.5 h' 1 Mpc. 
Transforming from Bj to H gives Mj Lh ~ 40 h M Q / L & for their data. This value is roughly 
half of our caustic estimate. Because our measured masses are almost exactly the same, the 
difference in M/L is due to their larger luminosity. Our luminosity determination is based 
on CCD photometry and agrees well with the optical CCD photometry of Kashikawa et al. 
(1995); Girardi et al.'s analysis is based on digitized and calibrated survey plates from the 
COSMOS/UKST Southern Sky Object Catalog and is discrepant with Kashikawa et al. 

Cluster mass-to-light determinations at H are scarce. We thus transform to R and V 
using the R — H colors above along with V — H ~ 3 for early-type galaxies and (V — H) Q = 
1.37. The mass-to-light ratios within 1.5 hr 1 Mpc are then equivalent to: 

M caustlc /L R = 374 ±59/i M Q /L Q 

M caustic /L v = 369 ± 59 h M Q /L Q . 
M virial /L R = 579 ± 119 h Mq/Lq 
M virial /L v = 572 ± 117 h Mq/Lq 

The caustic mass-to-light ratio is consistent with most previous measurements of cluster 
M/L ratios; typical values using virial masses are M/L R ~ 300 h M Q /L Q (Carlberg et al. 
1996; Girardi et al. 2000), with a range of 200 - 600 h M & /L &)R (Dressier 1978). Our virial 
M/L determination is consistent with the high end of this range. Estimates of M/L using 
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X-ray masses tend to be somewhat lower. Hradecky et al. (2000) find M/ Ly = 154 — 468 h 
solar units for seven nearby Abell clusters and one group, and ROSAT observations of eleven 
groups and clusters yield M/ L v = 200 — 300 h M Q /L Q (David, Jones, & Forman 1995). 

5. DISCUSSION AND SUMMARY 

We discuss a large, deep imaging survey of a galaxy cluster in the near-infrared, covering 
a 1.35 deg 2 (9 h~ 2 Mpc 2 ) region of A1644. These data allow us to determine the IR LF at 
H to roughly + 3, deeper than most previous infrared surveys of comparable size. We 
also acquired spectra for 155 galaxies (including ~ 141 cluster members) and use these to 
compute M/L H for the cluster. The results are: 

1. A grouping of cluster members by spectral type reveals evidence at the 96-97% con- 
fidence level for the segregation of emission and absorption galaxies by both velocity 
and core radius. Although only marginally statistically significant, these detections of 
segregation with respect to spectral type are in the same direction as the morphologi- 
cal segregation observed in other clusters (Adami, Biviano, & Mazure 1998; Melnick & 
Sargent 1977; Biviano et al. 1992). In contrast with an earlier analysis by Dressier & 
Shectman (1988b), there is no significant evidence for substructure on any scale. Our 
50% larger sample of cluster members is responsible for the difference. 

2. The velocity dispersion of the cluster varies little with radius for R < 1.5 hr l Mpc. 
Near the center, a is slightly larger than at the periphery, provided that we exclude 
three high- velocity galaxies at large clustercentric distance which may be background. 

3. The luminosity distribution is well-fit by a Schechter function with a = — 1.14±0.08 and 
Mft = —24.3 ± 0.2. These parameters agree with previous infrared LF measurements 
of clusters and field galaxies. In particular, an //-band determination of the Coma 
LF (Andreon & Pello 2000) yields an identical slope and characteristic magnitude to 
within the errors. The shape of the A1644 LF is consistent with several infrared cluster 
LFs (Trentham & Mobasher 1998; Barger et al. 1998; De Propris et al. 1999). However, 
due to small number statistics and the relatively flat slope of the LF at H ~ 13 — 13.5 
we are unable to rule out the presence of a dip similar to the ones observed in optical 
and infrared cluster LFs (Koranyi et al. 1998; Andreon & Pello 2000). 

4. We compute dynamical masses of A1644 based on the virial theorem and on the use 
of caustics. The caustic technique yields masses which have smaller errors, are more 
robust, and agree better with X-ray data. These masses, combined with the i/-band 
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photometry, yield an M/L H ~ 82 - 127 h M Q /L Q within R = 1.5 h' 1 Mpc. This 
mass-to-light ratio is one of the first determinations in the infrared, and corresponds 
to ~ 374 — 579 h M Q /L Q at visible wavelengths. That the infrared M/L agrees with 
previous optical determinations suggests that in some clusters biases by factors such as 
dust and short-lived star formation events are a small effect compared with the errors 
in the determination of mass-to-light ratios. 
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Fig. 1. — Magnitude difference versus dimmer magnitude for objects observed in two succes- 
sive dithered frames on 1997 April 17. At MAG_BEST ~ 16 the rms magnitude difference 
is 0.12 mag. 
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Fig. 2. — a) Velocity distribution of all 155 galaxies observed with FAST, b) Cluster sample 
of 141 galaxies in the range 10, 000 < cz < 20, 000 km s~ x . c) 123 cluster members for which 
i/-band photometry is available. In the bottom two plots, the dotted vertical line is the 
velocity of the cD while the solid line is the mean radial velocity of the cluster. 
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Fig. 3. — Sky positions of the 141 cluster members. Right ascension and declination are in 
decimal degrees. Open circles represent Abs galaxies and filled squares show Em systems. 
The cD, an absorption galaxy, is indicated by a cross. 
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Fig. 4. — Velocity distribution of the cluster sample separated by spectral type. Solid 
histogram: 123 Abs galaxies. Dotted histogram: 18 Em galaxies. Solid and dashed vertical 
lines indicate the mean velocity of the cluster and the velocity of the cD, respectively. 
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Fig. 5. — Top: Value of the Dressler-Shectman statistic A for varying subgroup size n. 
Bottom: Probability that A > A obscrvcd , based on 5000 Monte Carlo models for each value 
of n. With n — 11 our analysis does not confirm the detection of substructure by Dressier 
& Shectman (1988b). 
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Fig. 6. — Radial velocity profile of Em (top) and Abs (bottom) galaxies. The horizontal 
dotted line shows the mean velocity of the cluster. There are three circled galaxies which 
may be outliers. 
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Fig. 7. — Velocity dispersion for groups of 11 galaxies moving outward from the cD. Error 
bars show 68% confidence levels for uncorrected samples. Top: Includes all 141 galaxies 
with 10,000 km s -1 < cz < 20, 000 km s -1 . Bottom: Excludes the three high- velocity outliers 
circled in Figure 6. 
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Fig. 8. — Caustics and cluster members in position-redshift space. The caustic curves A(R) 
are shown as solid lines with la error bars. 
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Fig. 9. — Virial and caustic masses illustrating the inherent uncertainties in each technique. 
Solid lines: caustic mass calculated with smoothing parameter q=10, 25, and 50. Dashed 
line: virial mass including 141 galaxies with 10,000 km s -1 < cz < 20,000 km s _1 . Dotted 
line: virial mass using only galaxies which lie within the q=25 caustics. 
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Fig. 10. — Radial mass profile of A1644. Solid line represents the caustic estimate; dotted 
lines show upper and lower envelopes of the virial estimate with surface term correction. 
Crosses show the X-ray determinations from Ettori, Fabian, & White (1997). 
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Fig. 11. — Differential iJ-band luminosity function. Dotted histogram: all objects classified 
as galaxies by the SExtractor CLASS_STAR parameter. Dashed line: background counts 
modeled as a 10°- 67m power law. Filled squares with error bars: number of cluster members 
in each bin after background subtraction. Solid curve: best-fit Schechter function. 
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Fig. 12. — Mass-to-light ratio as a function of projected radius. Solid and dotted lines 
represent M/Lh computed with virial and caustic mass estimates, respectively. We show 
characteristic error bars at several locations along each curve. 
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Table 1. Properties of Photometric Sample 









cz 


RA (J2000) 


Dec. (J2000) 


rn H 


km s -1 


(1) 


(2) 


(3) 


(4) 


12 54 50.9 


-17 01 28.6 


14.46 




12 54 51.4 


-16 57 34.6 


15.19 




12 54 57.9 


-17 00 59.0 


15.15 




12 55 06.2 


-17 42 23.0 


13.97 


14971 



Note. — The complete version of this tabic is in 
the electronic edition of the Journal. The printed 
edition contains only a sample. 
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Table 2. Properties of Spectroscopic Sample 







cz 


O'cz 






DS 


RA (J2000) 


Dec. (J2000) 


(km s- 1 ) 


(km s- 1 ) 


Type 


m H 


Number 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


12 58 35.9 


-18 04 11.8 


12774 


33 


Abs 


13.18 


1 


12 56 10.2 


-17 59 44.7 


16712 


36 


Abs 


14.20 


2 


12 56 54.5 


-17 56 51.2 


13083 


16 


Em 


12.16 


6 


12 52 20.9 


-17 39 16.0 


13306 


34 


Abs 




9 


12 57 20.7 


-17 52 36.4 


15140 


33 


Abs 


13.86 


12 



Note. — The complete version of this table is in the electronic edition of the Journal. 
The printed edition contains only a sample. 



